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ABSTRACT 

iMatsuoka fe Kawaral (|2010f ) showed that the number density of the most massive galaxies 
(logM /Mq = 11.5 — 12.0) increases faster than that of the next massive group (log M/Mq = 
11.0 — 11.5) during < z < 1. This appears to be in contradiction to the apparent "downsiz- 
ing effect". We attempt to understand the two observational findings in the context of the hier- 
ar chical merger paradigm usi ng semi-analytic techniques. Our models closely reproduce the result 
of lMatsuoka fc Kawaral (|2010f ). Downsizing can also be understood as larger galaxies have, on average, 
smaller assembly ages but larger stellar ages. Our fiducial models further reveal details of the history 
of the stellar mass growth of massive galaxies. The most massive galaxies (log M/Mq — 11.5 — 12.0 
at z=0), which arc mostly brightest cluster galaxies, obtain roughly 70% of their stellar components 
via merger accretion. The role of merger accretion monotonically declines with galaxy mass: 40% for 
log M/Mq = 11.0 - 11.5 and 20% for log M/Mq = 10.5 - 11.0 at z = 0. The specific accreted stellar 
mass rates via galaxy mergers decline very slowly during the whole rcdshift range, while specific star 
formation rates sharply decrease with time. In the case of the most massive galaxies, merger accretion 
becomes the most important channel for the stellar mass growth at z ~ 2. On the other hand, in-situ 
star formation is always the dominant channel in galaxies. 

Subject headings: galaxies: evolution - galaxies: elliptical and lenticular, cD - galaxies: formation 
galaxies: stellar content 



1. INTRODUCTION 

Dynamical realizatio ns based on the c oncordance 
ACDM cosmology (e.g. iSpergel et al.l 120071 ) have been 
remarkably successful at r eproducing large-sca le struc- 
tures in the universe (e.g. iSpringel et all 120061) . In the 
common perception of this paradigm, large galaxies grow 
hierarchically through numerous galaxy mergers that fol- 
low mergers between dark halos. 

Supporting this view, galaxies with dist urbed fea- 
tures have frequently be en witnessed (e.g. lArpI 119661 : 
ISchweizer fc Sc itzcr 1988). With the advent of deep and 
wide field surveys, the arguments could be investigated in 
greater detail. Recent studies based on ultra-deep imag- 
ing data are particularly not eworthy. From the u r = 28 
mag arcsec -2 deep images, Ivan Dokkuml (p005T ) found 
that about 50% of red bulge-dominant galaxies in field 
environments show tidal debris. Merger galaxy fraction 
has been fo und to be almost as high in cluste r envi - 
ron ments fSheen et al.ll2012f ). iKavirai et all (|2007f ) 
and lKavirail(|2010H claimed that residual star formation 
found in a large fraction (~ 30%) of local massive early- 
type galaxies is related to mergers or interactions. As 
observation techniques reach deeper and hidden nature 
of galaxies, "peculiar" is no longer a synonym of "rare" . 

High rcdshift surveys allow us direct investigation 
of the role of galaxy mergers and interactions on 
galax y evolution. Using the K-b and Hubble dia- 
gram, [Ar^'on-SalamancjieLal. (1998) found that bright- 
est cluster galaxies(BCGs) have increased their mass by 
a factor of two to four with no or negative luminos- 
ity ev olution durin g < z < 1. From the GOODS 
fields, iBundv et all (|2009h found that the pair fraction 
of massive (logM/M Q > 11.0) red spheroidal galaxies is 
higher than that of less massive systems (logM/M Q ~ 



10.0). In addition, the pair fraction unaccompanied by 
star formation increases with time. They concluded 
that massive galaxies grow primarily through dry or 
minor mergers, at least at z < 1. Almost simultane- 
ously, fr om the UKIDSS and the S PSS II Supernova 
Survey, IMatsuoka fc Kawaral (|2010t hereafter MK10) 
found that, during z < 1 the number density of the 
most massive galaxies (log M/Mq = 11.5 — 12.0) in- 
creased more rapidly than that of the next massive group 
(log M/Mq = 11.0 - 11.5). They also showed that more 
massive systems have a lower blue galaxy fraction than 
less massive systems, and that the fraction decreases with 
time. All these observations seem to imply that, during 
z < 1, massive galaxies are mainly brought up via merg- 
ers. 

Not all observations naively support the hierarchi- 
cal merger picture. The tight color-magnitude rela- 
tion found among early-type galaxies is more simply, al- 
beit not exclus ively, explained by monolithic formation 
scena rios (e.g. I Bower et al.l Il992t iKodama fc Arimotol 
|1997[ ). Furthermore, the "downsizing" effect, where 
larger galaxies appear to be older and thus suspected 
to have formed ea rlier, seems to be inconsistent with 
the new paradigm (ICowie et alj|l996t IGlazebrook et ahl 
12004 ICimatti et all 120041) . On the face of it, the in- 
consistency seems a counter-evidence of the hierarchi- 
cal picture; however, some studies have pointed out 
that downsizing can be understood by the hierarchical 
paradig m reasonably e n ough . Based on semi-analytic 
models, IDe Lucia et al.l ([2006D showed that the star for- 
mation rates of all the progenitors of more massive galax- 
ies peak earlier and decrease faster than those of less 
massive galaxies. In fact, downsizing could be a nat- 
ural n3Su]t_of_the hierar chical clusterings of dark ha- 
los. INeistein e t al.l (|2006f ) demonstrated that downsiz- 
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ing appears in the total (combined) mass evolution of 
all the progenitors of a dark matter halo. If the stel- 
lar mass of a galaxy correlates w ith the mass of its dark 
matter halo (jMoster et alJ[201Cf ). downsizing may very 
well originate from the bottom-up assembly of dark mat- 
ter halos. Using semi-analytic approaches, many stud- 
ies have shown that the assembly and formation his- 
tory of stellar components, especially in massive, red 
and luminous, or elliptic al galaxies, can be different in 
the h i erarchical Unive r se (iKauffmann 1996; B augh et alJ 
19961: iDe Lucia et all 120061 : iDe Lucia fc Blaizotl 120071 : 
Almeida et al.ll2008D . More direct hydrodynamic simu- 
lations perfor med for smaller volumes have yielded con- 
sistent results ()Oser et al.ll2010t lLackner et al.ll2012ft . 

In this study, we look further into the details of the as- 
sembly history of stellar components in massive galaxies. 
We use semi-analytic approaches because these are more 
effective for constructing models of a large volume of the 
universe. Motivated mainly by the empirical results of 
MK10, we investigate the factors that drive the growth 
of galaxy stellar mass as a function of time and mass. 



2. MODEL 

We have developed our own semi-analytic model for 
galaxy formation and evolution. In this section, we 
briefly introduce physical ingredients, with more focus on 
the prescriptions that play important roles in this study. 

2.1. Dark Matter Halo Merger Trees and Galaxy 
Mergers 

In the context of the two-step galaxy formation the- 
ory (j White fc Reesl[T978| ) , the first step involves the con- 
struction of dark halo merger trees. We ran dark mat- 
ter N -body volume simulations using the GADGET-2 
code (|Sprin gcl 2005|). The simulation was performed us- 
ing the standard ACDM cosmology parameters derived 
from the WMAP 7-year observations, Vl m = 0.266, JIa = 
0.734, cr s = 0.801, »h = .0449, and H = 71 km s" 1 
Mpc -1 (jJarosik et al.l [20TTI ) . The periodic cube size of 
the simulation was 70/i _1 Mpc on the side with 512 3 colli- 
sionless particles, mass of a particle was 1.9x 1O 8 /i _1 M0, 
and mass resolution of a halo was ~ 10 10 h~ 1 M ( 7 ) . We 
identify halo structures in the si mulation box , usin g 
a halo finding code developed by iTweed et al.l (2009 ). 
which is based on the AdaptaHOP (jAubert et al.ll2004l ). 
Halo merger trees were generated by backtracing the in- 
fall histories of subhalos with increasing rcdshift. We 
calculate galaxy merger timescales, using the positional 
information of subhalos extracted from our N-body sim- 
ulation. In this study, we assume that the mass distribu- 
tion of dark matter halos ba sically follows the N avarro- 
Frank- White (NFW) profile (jNavarro et al.lll996[ ). 

Subhalos can disappear before arriving at central re- 
gions if they are heavily embedded in their host halo den- 
sity profiles or their mass decreases below the resolution 
limit of our numerical simulation. We take these numeri- 
cal artifacts into consideration. Thus, once a halo enters 
into a more massive halo, we additionally calculate its 



merger timescale, t^ 



mula suggested by Uiang et al 



using the fo llowing fitting for- 
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ln[l + (Mhoat/Msat)] Vc 
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where e is the circularity of the orbit of a satellite halo, 
C is a constant, approximately equal to 0.43, Mh os t is 
the mass of a host halo, M sat is the mass of a satellite 
halo, i? v ir is the virial radius of a host halo, and V c is 
the circular velocity of a host halo at i? v ir- If a subhalo 
disappears within 0.1i? v ; r of its host halo, a galaxy in 
the subhalo is regarded as having merged with its cen- 
tral galaxy. On the other hand, if it disappears out- 
side 0.1i? v ir, we assume that the galaxy in the subhalo 
merges with its central galaxy at i mC rgc given by Eq. [1] 
after the subhalo becomes a satellite of its host halo. 
In that case, we consider dynamical friction to analyt- 
ically compute the positions and velocities of subhalos. 
W e adopt the dynamical frict ion prescription introduced 
bv lBinnev fc Tremaind (pOOl : 




v 



cxp 



(2) 



where M sat is the mass of a subhalo, which is initially 
defined as the mass at a previous time step after which 
the subhalo cannot be resolved in the N-body simulation 
anymore, r is the distance between the subhalo and the 
center of its host halo, In A is a C oulomb logarithm wit h 
A = 1 + M host /M 8at adopted by ISpringel et al.1 (|2001D , 
V c is the circular velocity of the host halo at the virial 
radius, and v is the orbital velocity of the subhalo. 

During the orbital motion of a satellite halo in its host 
halo, the dark matter of the satellite halo is stripped due 
to dynamical friction. If the satellite is resolved in an 
N-body simulation, dark matter of the satellite would be 
naturally stripped. On the other hand, if the halo is not 
resolved in the simulation but considered to prbit around 
its host halo, stripping due to dynamical friction should 
be computed analytically. We evaluate the amount of 
dark matter stripped by dynamical friction by adopting 
the concept of the sphere of influence (r so i) within which 
dark matter particles are b ound to the satellite halos, 
using the following formula (]Battin|[i987l ): 



A/host (< r) 



l 1 + 6 cos 

(1 + 3 cos' 8) + 0.4 cos e 

'1+3 cos 2 
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where r is the distance between the centers of the satel- 
lite and its host halos, M satj tot is the total (baryon+dark 
matter) mass of the satellite halo, Mhost (< r ) is the to- 
tal mass of the central halo within r, and 6 is the angle 
between the line connecting the particle to the center of 
the satellite halo and the line connecting the centers of 
the satellite and the host halos. During a time step, St, 
we assume that a satellite halo loses SM sat = M sat (r > 
r S oi)o~t/td yni where £d yn is the dynamical timescale of the 
satellite halo, and M sat (r > r so j) is the mass of dark 
matter outside r SO [ of the satellite halo. 

We assume that stellar components in satellitcQ galax- 
ies merging with their hosts constitute the bulge corn- 
In this paper, galaxies that are not the central one in a halo are 
all "satellite" . Only one galaxy is qualified as the central galaxy of 
a halo and all the rest, regardless of brightness, arc satellites. 
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ponent of the host. If the mass ratio of bary- 
onic mass (m co id + m*) between merging galaxies, 
ra SCC ondaiy/«Vimary, is greater than 0.25, then it is as- 
sumed that all the stellar components of the host galaxy 
quickly become bulge components of the remnant, as 
well. 

Empirical studies have shown that intra-cluster light 
originates from the extended diffuse stellar components 
of the brightest galaxies in groups or clusters (e.g. 
Feldmeier ct al. 200l iGonzaleTc t al. 20051 iZibetti et al.l 



20051 ). They suggested that diffuse stellar components 
are the stars scattered from satellite galaxies during 
tidal stripping or mergers into central galaxies. It has 
been suggested that 10-40% of stellar components in 
satellite galaxies tur n into diffuse stellar components in 
each galaxy merger (jMurante et al.1 120041: iMonaco et"al~l 
I2006T ) . We adopt the value of 40% in this study because it 
resulted in the best reproduction of empirical data. The 
amount of stellar mass that a central galaxy acquires via 
merger is (l-/ S catter)Af*, 8 at, where /scatter is the fraction 
of scattered stellar components and -M* jSat is the stellar 
mass of a satellite galaxy. 

2.2. Gas Cooling, Star Formation, and Recycling 

We assume that the baryonic fraction in accreted 
dark matter follows the global baryonic fraction, f2b/^ m - 
Baryons are accreted onto dark halos and shock-heated 
to become hot gas components. 

Gas accretion onto a galactic disk plane via atomic 
cooling of hot gas is calculated based on the model pro- 
posed by iWhite fc Frenkl (|1991[ ). The cooling timescale 
at distance r from the center of a halo is estimated by 
the following formula: 



cool 



(r) = 7T 



p g {r)kT 



2 p(Z, T)mpn e niA(Z , T) ' 



(4) 



where p g is the gas density within radius r of a dark 
matter halo, T is the temperature of gas, assumed to 
be the virial temperature in the model, Z is metallicity, 
p(Z, T) is the mean molecular weight of gas with Z and 
T, mp is the mass of a proton, n e is the number density of 
electrons, is the number density of ions, and A(Z,T) 
is the cooling function with Z and T. The values of 
p( Z,T), n R , n,;, and A(Z , T) ar e determined by referring 
to ISutherland fc Dopital (|1993l ). We do not include self- 
consistent chemical evolution in our calculation. Instead, 
we take the metallicity of hot gas components in halos as 
a constant, 0.3Zq, which is comparable to the metallicity 
of ob served clusters with various masses ([Arnaud et al.l 
I1992H . Because the cooling function, A(Z,T), is sensi- 
tive to metallicity, our results should not be taken too 
literally. However, relative analysis, a main tool of this 
study, is affected little by the details in the treatment of 
chemical evolution. 

It is assumed that the gas density follows a sin- 
gular isothermal profile truncated at i? V ir^ Pg( r ) = 



V 



Substituting the formula for p g (r) 



in Eq. |4] and adopting the dynamical timescale of a 
halo, tdym as t coo i, one can derive the cooling radius, 
fcooi, within which hot gas can cool within t coo \. For 
the case of r coo i > R v i T , the cooling rate is rather re- 
strained by the free-fall rate than the cooling rate, so 
that m cool = m ho t/(2£ C ooi)- In contrast, if r cool < i? vir , 



^cool — ^hot^cool/(2^vir^cool ) • 

In our model, stars can be formed through a quiescent 
mode, in which cold gas turns into disk stellar compo- 
nents via gas contraction on a disk, or a burst mode, 
which is induced by galaxy mergers. Star formation rate 
in the quies cent mode are de l ineate d by a simple law 
proposed bvlKauffma nn et al.l (|1993|) as follows: 



TO* 



TO co ld 

1 t, — 

^dyn^gal 



(5) 



where a is the empirically-determined star formation ef- 
ficiency, TO co i(j is the amount of cold gas, and idyn.gai is 
the dynamical timescale of cold gas disk assumed to be 

O.lidyn- 

Ob servations (e.g. iBorne et al.1 120001: iWoods fc Gelled 
I2007T) a nd hydrodynam ic simulations of galaxy merg- 
ers (e.g. ICox et al~ll2008l hereafter C08) have shown that 
galaxy mergers can give rise to rapid star formation. We 
follow the conventional treatment: stars formed in the 
quiescent mode belong to a galactic disk, while stars born 
in the burst mode become bulge components. We adopt 
th e prescription for merge r-induced starbursts described 
in ISomerville et all ([20081 . hereafter S08), which formu- 
lates the prescription based on C08. S08 defines burst 
efficiency, eburst, to parameterize the fraction of the cold 
gas reservoir involved in a merger induced starburst as 
follows: 



eburst — eburst, 0/^ 



(6) 



where p is the mass ratio between a host galaxy and 
its merger counterpart, 7burst is the bulge-to-total mass 
ratio(B/T) of the host galaxy, and eburst, o is the burst 
efficiency fitted by the following formula: 

eburst.o = 0.60[Kh-/ (kms- 1 )] - 07 ^ + ffiDOsr - 17 

(l + /g) 0.07 (1 + z) 0.04 ; (7) 

where Vv U - is virial velocity, (jEOS is the effective equation 
of state of gas, / g = TO co id/(TO co id + to*) is the fraction 
of cold gas, and z is the redshift when the disks of pro- 
genitor galaxies are constructed, qeos was suggested to 
parameterize the multiphase nature of ISM: <7eos = in- 
dicates an isothermal state, and qeos = 1 represents the 
fully pressurized multiphase ISM. In this study, we adopt 
9eos = 1 at which gas is dynamically stable, so that star- 
bursts are more suppressed than the case of (?eos = 0. 
The redshift dependency of eburst, o is very weak, and thus 
we assume (1 + z) 04 ~ 1. The burst timescale, Tburst, is 
also formulated by S08 as follows: 

Tburst = 191Gyr[Kh7 (kms" 1 )]- 1 - 8 ^ + ^eos) 2 ' 58 

(i + / g r a74 (i + ^r - 16 - (8) 

The mass ratio, p, is the ratio of the total mass of central 
regions (mDM,core+TO^+7n co id) of a host galaxy to that of 
its merger counterpart. Following S08, we calculate the 
core mass of a dark matter halo, todm.coio = todm(?" < 
2r s ), where r s = i? v ir/cNFW- The concentration index of 
the Navarro-Frenk- White profile, c nfWi is derived base d 
on the fitting function suggested bv lMaccio" "eTal](|2007h . 
The parameter 7b U rst is determined by the B/T of a 
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host galaxy as follows: 



Tburst 




B/T < 0.085 
0.085 < B/T< 0.25 
0.25 < B/T 



(9) 



C08 showed that the burst efficiency of a host galaxy 
with a high B/T is lower than that of a galaxy with 
a lower B/T, because more massive bulges stabilize the 
galaxies and reduce the burst efficiency more effectively. 
Because C08 demonstrated that mergers with mass ra- 
tios below 1:10 are not associated with starbursts, we 
assume eburst = if fx < 0.1. With the ingredients, 
the amount of stars born in burst modes is calculated as 
tburst = eburst TOcoid- We assume that mburst turns into 
stars for r burs t in a uniform rate, m bu rst = "i bu rst /n>urst- 
While merger-induced starbursts occur in a galaxy, the 
quiescent mode also still goes on in our models. 

In our model, the recycling of stellar mass loss is con- 
sidered in great detail. We compute the mass loss of 
every single population at each epoch after a new stel- 
lar population is born. The mass loss of a single pop- 
ulation is calculated as follows. We adopt the Scalo 
initial mass function ([Scalol 119861) . The lifetime of a 
star with mass M, tm, is computed by a broken-power 
law dFerreras fe Silkl I2000D , which is obtained f rom the 
data of iTinslevi ([19801) and lSchaller etail (fl992T) . 



T M (Gyr) 



9.694(M/M )- 2 - 762 
O.O95(M/M )-°- 764 



M < 1OM 
M > 10M, 



■ 



(10) 



At the end of its lifetime (after tm elapses), a star returns 
most of its mass into space, leaving small remnants such 
as a white dwarf, a neutron star or a black hole. The 
remnant mass of a star with mass M, ojm, is suggested 
by Ferreras & Silk (2000) as follows: 



•0.1(M/M Q )4 
1.5 

M © io.61(M/M ) 



0.45 



M < 1OM 
lOAf < M < 25Me(ll) 
13.75 M > 25M . 



In this study, we simply assume that half of the mass loss 
returns to cold gas components and the rest becomes hot 
gas. 

2.3. Environmental Effect 

The Chandra X-ray Observatory revealed that massive 
satellit e galaxies in nearby clusters have hot gas compo- 
nents (jSun et al.l 120071 : Ueltema et al.l I2007D , in contra- 
diction to expectations based on an instantaneous hot 
gas stripping sc enario for satellite galaxies in cluster en- 
viron ments (e.g. lKauffmann et al.lll99"9l : lSomerville et al.l 
l2008f ) . The old assumption predicted a higher fraction of 
passive satellites in large halos than observed, known a s 
the satellite o verqu enching problem (jKimm et^[2009h . 
IKimm et al.l (|2011D showed that a gradual, rather than 
instant, and more realistic stripping of the hot gas reser- 
voir relieves the above-mentioned problem to some de- 
gree. Therefore, we implemented the gradual hot gas 
stripping of satellit e galaxies in our mo del by considering 
tidal stripping (see Kimm e t al.ll20l"T[ ) and ram pressure 
stripp ing that uses the prescriptions of IMcCarthv et al~l 
(|2008f ) . which had been modified for semi-analytic mod- 
els by t^^eFil] (HOOl) . However, we are still missing 
a realistic prescription on cold gas strip ping. Empiri- 
cal evidence for cold gas stripping is clear (|Vollmer et al.l 



120081: iChung et al.|[2008allH ). and a theoretical study us- 
ing a sem i-analytic approac h shows its effect in galaxy 
evolution dTecce et al.ll2010l ). Thus, it should be consid- 
ered in our model in due course. 



2.4. Feedback Processes 

Feedback mechanisms have been introduced into 
galaxy formation theory to reconcile the discrepancy be- 
tween galaxy and dark matter halo mass functions. Hy- 
drogen atoms, neutralized at the recombination, may be 
rcionized at later epochs (z < 10) due to background 
high-energy photons. Thi s mechanism may suppress the 
growth of small gala xies (|Gnedinll2000t ISomerville! [20021: 
IBenson et al.1 12002alrb[ ) . Supernova feedback is thought 
to be effective at dist urbing the growth of s mall galaxies 
by ej ecting cold gas (|White fe Reesl 119781: iDekel fe Silkl 
I1986I) . and AGN feedback is considered m ore effective 
in massive galaxies with a large bl ack hole (|Silk fe Reesl 
[llMIScliawinski et al j[2fM [2001 . 

We u tilize the reionization prescription of lBenson et al.l 
(|2002bD . The prescription allows an inflow of baryons 
into a dark halo via accretion of dark matter when 
V-vix > K-oionization throughout the age of the Universe, 
where ^ionization is the suppression velocity of reioniza- 
tion . If a halo has a lower value of Kdr than the criterion, 
the inflow of baryons is allowed only at z > ^reionization, 
where ^reionization is the suppression redshift of reioniza- 
tion. In this Study, we adop t Vreionization = 30 km s _1 

and ^reionization = 8 following IBenson et al.l (|2002bD . 

We follow the prescriptions of S08 for supernova feed- 
back. These prescriptions take into account not only the 
amount of reheated gas but also the fraction of reheated 
gas blown away from halos. The reheating rate of cold 
gas due to supernova feedback is formulated as follows: 



TO rh 



J3N 



/ 200km f 



— 1\ a rh 



isk 



(12) 



where 6q and a rn are free parameters, Vdisk is the ro- 
tational velocity of a disk, and to* is the star formation 
rate. S08 assumes that the rotational velocity of a disk is 
the same as the maximum rotational velocity of the DM 
halo. We calculate the fraction of the reheated gas that 
has enough kinetic energy to escape from the halo as 



/ eject 



1.0 



v. 



eject 



(13) 



where a e joct = 6 and V e j e ct ~ 100 — 150km s~ x . In the 
case of a satellite halo, a fraction of reheated gas by su- 
pernova feedback, /ejectWrh, is ejected from the satellite 
and added to the hot gas reservoir of its host halo. On 
the other hand, it is assumed that the ejected gas from 
main halos is diffused through inter-cluster medium. 

We take the quasar-mode and radio-mode AGN feed- 
back into account in our model. It has been sug- 
gested that the quasar mode is induced by an in- 
flow of cold gas into the central super-massive black 
hole (SMBH) of the central galax y during major merg- 
ers (jKauffmann fc Haehneltl [20001 ). The increasing mass 
of the SMBH via the accretion of cold gas can be ex- 
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Fig. 1. — The cosmic star formation history. The gray dots 
with error bars indicate the empirical cosmic star formation his- 
tory l lPanter et al.l [20071) . The solid line shows the cosmic star 
formation history derived from the fiducial model. The dotted line 
displays the contribution to the cosmic star formation history from 
merger-induced starbursts. 



pressed as 



Am B H,Q = 



/]3H m cold 



1 + (280km S"VKi 



(14) 



where / BH is the efficiency of gas accretion. In 
this study, we take the modified parameter proposed 
bv lCroton et all (j2006D : / BH = / B H(A/ sat /M host ), where 
/b h is the original form o f the p arameter introduced 
bv iKauffmann fc Haehneltl <|2000T >. M host is the mass 
of a host galaxy, and M sat is the mass of the host 
galaxy's merger counterpart. In this study, we assume 
th at the lifetime of the quasa r m ode, toso, is 0.2Gyr 
as IMartini fc Weinberg] (|2001l ) and IMartini fc Schneider! 
(|2003f ) suggested tnm < 0.3Gyr. Thus, we have 
'tibh.q = Atobh,qAqso- It is generally thought that 
the quasar mode is caused by a high accretion rate of 
cold gas, resulting in a rapid growth of an SMBH. 

The radio-mode feedback releases low-Eddington-ratio 
energy through the accretion of hot gas distributed 
throughout the halo. Although the energy released from 
the radio-mode AGN is far less than that from the quasar 
mode, it is regarded that the radio mode supplies enough 
energy to the surrounding medium to interrupt gas cool- 
ing or to blow away (some of the) cold gas. We imple- 
ment radio-mo de feedback into ou r model following the 
prescription of ICroton et al.l (|2006[) : 



«^BH,R — « AGN 



U, 



/ ra-BH \ / /hot \ / 

V10 8 M Q y V 0.1 / \ 200km s" 



(15) 



where kagn is a free parameter with units of Mq yr^ 1 , 
mBH is the black hole mass, /hot is the fraction of hot 
gas with respect to the total halo mass, and V v „ is the 
virial velocity. 

We assume that the amount of energy generated by 
the accretion of gas into the SMBH is given as follows: 



L B h = rjmBHC 2 = ^(ttibh.q + "ibh,r)c 2 



(16) 



where 77 = 0.1 is the standard efficiency of the conversion 
of rest mass to radiation, and c is the speed of light. The 
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Fig. 2. — The galaxy s t ellar mass functions in the local Universe 
derived bv lPanter eTaTl pOOTT ) from SDSS DR3 (gray shade) and 
the fiducial model at z=0 (black solid line). The thickness of the 
empirical data indicates the error range. 



reduced cooling rate of gas is computed by 



TO cool — m cool 

where the minimum of m j 



o.wx 

is set to be zero. 



(17) 



2.5. Model calibrations 

Our models are based on the conventional techniques 
and ingredients used in up-to-date semi-analytic mod- 
els; hence, the output is not particularly noteworthy 
compared to other successful models. Our models 
roughly match the global star formation history, galaxy 
mass functions, black hole mass versus bulge mass re- 
lation, etc. Figure [1] and [5] display comparisons of 
the cosmic star formation history and the galaxy stel- 
lar m ass functions in th e local Universe from empirical 
data ijPanter et al.1 120071 ) and our fiducial model. While 
there still is a large room for improvement, we decide to 
focus on the mass growth histories of massive galaxies. 



3. EVOLUTION OF GALAXY NUMBER DENSITY 

MK10 presented a rapid growth of massive galaxies 
since z = 1, using the United Kingdom Infrared Tele- 
scope(UKIRT) Infrared Deep Sky Survey (UKIDSS) and 
the Sloan Digital Sky Survey (SDSS) II Supernova Sur- 
vey. Figure [3] shows the number density evolution of the 
most massive (logM/M© = 11.5 — 12.0) and the next 
massive (log M/M® = 11.0— 11.5) galaxies in the empir- 
ical data derived by MK10 and from our fiducial model at 
each redshift. The empirical data clearly show that the 
number of the most massive galaxies rapidly increases 
between z = 1 and 0, while the next massive group expe- 
riences a milder evolution. The reproduction of the data 
by our fiducial model looks reasonably good. We also 
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Fig. 3. — Number density evolution of massive galaxies as a 
function of the age of the Universe. The red represents the evo- 
lution of the most massive galaxies (log M/Mq = 11.5 — 12.0), 
the blue indicates that of the next massive group (log M/Mq = 
11.0—11.5), and the green indicates that of the third massive group 
(log M/Mq = 10.5 — 11.0) at each redshift. The diamonds with er- 
ror bars come from empirical data in M K10 and the square s with 
error bars are measurements taken from ICole et al.l H200JJ ) . The 
solid lines present the predictions of the semi-analytic model. 

present the number density evolution of the third mas- 
sive group (log M/Mq = 10.5 — 11.0), which makes our 
models "super-L* galaxies". We define "relative num- 
ber density growth rate", T, as the ratio of the speeds in 
number density evolution of the two most massive groups 
of galaxies as a function of observational limit in redshift, 
as follows: 

„ n most {z = 0)/n most (z) 
n nC xt(z = 0)/n ncxt (z) 

where n most (z = 0) and n nex t(z = 0) are the number den- 
sities of the most massive and next massive groups at 
z = 0, and n most (z) and rt noxt (z) are the number den- 
sities of the two groups at a redshift, respectively. For 
example, MK10 compared the number density evolution 
of the two mass groups between redshift and 1 , in which 
case the observational limit is 1 and the relative number 
density growth rate becomes 3. We present their observa- 
tions and our models in Figure [4] The MK10 data point 
should be compared with our model for the most and 
the next massive galaxy groups (solid line). The model 
that compares the next massive group of galaxies with 
L* galaxies (dashed line) exhibits a similar but milder 
trend. Observational constraints are still weak but are 
at least roughly reproduced by the models. The relative 
growth rate is always greater than 1 in the models, which 
indicates that the number density of more massive galax- 
ies undergoes a faster evolution than that of less massive 
groups in the super-L* range. The number evolution is 
most dramatic when a comparison is made against the 
most massive galaxies. This is caused by the fact that 
the mass bin for the "most massive" galaxies has only 
influx from less massive galaxies, whereas the mass bin 
of less massive galaxies can have outflux to more mas- 
sive galaxy bins as well as influx from even less massive 
galaxy bins. 

Primarily motivated by MK10, we focus on the mass 
growth histories of super-L* galaxies. We divide the 
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Fig. 4. — Relative galaxy number density growth rate, T, be- 
tween z=0 and observational redshift limits. The solid line shows 
r of the most/next and the dashed line indicates that of the 
next/third in our model. The cross is derived from the empiri- 
cal data in Figure [3] 

model galaxies into three groups according to z = mass: 
Rank 1: log M/M = 11.5-12.0, Rank 2: log M/M Q = 
11.0-11.5, and Rank 3: log M/M© = 10.5 - 11.0, where 
Rank 3 roughly represents L* galaxies. From our simula- 
tion volume, we found 49 galaxies in Rank 1, 472 galaxies 
in Rank 2, and 2,188 galaxies in Rank 3. 



4. EVOLUTION OF MASSIVE GALAXIES IN 
MODELS 

4.1. Evolution of Stellar Mass in Galaxies 

In the hierarchical paradigm, a galaxy can have more 
than one progenitor. Progenitors of a galaxy can be di- 
vided into "direct" and "collateral" progenitors. A di- 
rect progenitor is the galaxy in the largest halo when 
a merger between halos takes place. While there can 
be numerous progenitors, there is only one direct pro- 
genitor at each epoch. Collateral progenitors are all the 
other galaxies that contribute to the final galaxy. In this 
concept, to build the evolutionary history of a galaxy, 
one should consider not only direct progenitors, but also 
merger counterparts or collateral progenitors. Figure [5] 
shows the average mass evolution of the direct progeni- 
tors (solid lines) and all (direct and collateral combined) 
the progenitors (dotted lines) of Rank 1,2, and 3 galax- 
ies. The mean stellar masses of the three groups at z = 
are 4.97 x 10 n M o , 1.55 x 1O U M , and 5.25 x 1O 1O M , 
respectively. Most of the Rank 1 galaxies in our volume 
are brightest cluster galaxies. 

As all the progenitors merge with each other, the dot- 
ted lines and the solid lines finally meet at z = 0. Stel- 
lar mass loss and the scattering of stellar components in 
satellite galaxies into diffuse stellar components, which 
takes place when mergers occur, lead to a gradual de- 
crease in the total mass during the evolution. For exam- 
ple, one can see a slight decline in the total stellar mass 
(red dotted line at the top) after z ~ 0.5. This effect is 
not clearly visible in Rank 2 and 3 galaxies in which star 
formation is more extended and mergers are less frequent 
than in Rank 1 galaxies. 

It is useful to have a definition of the formation red- 
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Fig. 5. — Average mass evolution of galaxies. The red, blue, and 
green represent three groups of galaxies: log M/Mq = 11.5 — 12.0, 
log M/Mq = 11.0 - 11.5, and log M/M = 10.5 - 11.0 at z = 0, 
respectively. The solid lines indicate the mean mass evolution of 
direct progenitors and the dotted lines show the mean mass evo- 
lution of all (dircct+collatcral) progenitors. The black horizontal 
dashed lines denote half the mean galaxy mass at z = of the 
groups. Ai, Aa, and A3, and the arrows indicate the epochs when 
the total stellar masses of all progenitors reach half of the final 
mass. Di, D2, and D3 with arrows show the epochs at which the 
direct progenitors of the three groups acquire half of their final 
mass. 

shift, Zf. We define it as the redshift at which half of 
the stellar mass at z = has been assembled. In the 
case of Rank 1 galaxies, half of the final mass is achieved 
at z^d ~ 0.7(denoted as Di in Figure [5]) in direct pro- 
genitors and at z^a ~ 3.0(Ai) when all progenitors are 
combined. Our models suggest (zf,D, Zf,A) = (0.9, 2.1) 
for Rank 2, and (zf,D, Zf,A) = (1-1 3 1-6) for Rank 3. 
Models exhibit a monotonic mass dependence of Zf,D 
and Zf^A in the sense that, with mass, z^d decreases 
while Zf : A increases. In other words, the mass of the 
direct progenitors of a more massive group grows more 
slowly, while its total mass of all the progenitors is as- 
sembled earlier than that of a less massive group. The 
evolutionary histories of direct progenito rs are opposite 
to the pat tern of cosmic do wnsizing. As iNeistein et al.l 
(|2006ft and lOser et"al] (|2010l) pointed out, however, if the 
growth histories of collateral progenitors of galaxies are 
also considered, downsizing would be a natural outcome 
of the hierarchical concept of galaxy formation. 

The difference in the growth history between the three 
groups can be understood in depth through Figure 
which presents the evolutionary histories of star forma- 
tion rates (SFRs). The mean SFRs contain the star for- 
mation histories of both direct and collateral progenitors. 
We show the best fitting log-normal function to the three 
SFR curves. The star formation rates of more massive 
galaxies peak earlier, as marked by SI, S2, and S3 in 
the figure, and decreases faster than those of less mas- 
sive grou ps. The same features were noted in an earlier 
study bv lDe Lucia et all (12006D . 

Figure [H] also show the redshift at which half of the to- 
tal cumulative star formation has occurred in the three 
groups of galaxies: (CI, C2, C3) = (2.0, 1.5, 1.4) in z. 
The general trend shown by these three values agrees 
with that of z^a discussed above. In a sense, it is these 
values, rather than z^a, that are closer to the general 
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Fig. 6. — Mean star formation histories of all the progenitor 
galaxies of Rank l(red), 2 (blue), and 3 (green). The gray dashed 
lines denote log-normal fitting. Si, S2, and S3 with arrows indicate 
the epochs when SFRs peak. Ci, C2, and C3 with arrows show 
the epochs at which the cumulative stellar mass reach half of the 
total stellar mass born by z = 0. 

definition of formation redshift. Again, the more massive 
galaxies are, the earlier they form their stars, consistent 
with downsizing. In conclusion, the observational finding 
of downsizing is a result of the hierarchical galaxy for- 
mation process, where more massive galaxies have larger 
stellar ages and smaller assembly ages (se e also e.g. 
iDe Lucia fc Blaizotll2007t iKavirai et al.ll2009f ). Both the 
larger stellar ages and the smaller assembly ages can be 
understood as a result of large-scale effect; that is, in 
a deeper potential well, progenitor galaxies and their 
stars form earlier, and many more galaxies participate 
in galaxy mergers for a long period of time. 

4.2. Origin of Stellar Components 

In this section, we investigate how stellar components 
are assembled into massive galaxies. Stellar components 
in a galaxy originate either from in-situ star formation or 
from "merger accretion" . Two modes of star formation 
are considered: "quiescent" mode and merger-induced 
"starburst". Stars in a galaxy can therefore have four 
different origins: (1) in-situ quiescent star formation, (2) 
in-situ starburst, (3) merger accretion of stars formed in 
quiescent mode, and (4) merger accretion of stars formed 
in burst mode. 

Figure [7] shows the decomposition of the four channels 
as a function of time for the three different mass groups. 
The top and bottom rows show the absolute stellar mass 
evolution and relative mass fraction evolution, respec- 
tively. 

Several remarkable features arc visible, more easily 
in the bottom rows. First, in-situ quiescent star for- 
mation is an important channel for the stellar com- 
ponents in massive galaxies. Its fractional contribu- 
tion is (30, 60, 80)% in Rank (1, 2, 3) galaxies, re- 
spectively. In-situ quiescent star formation takes place 
on a galactic disk, and thus, one may wonder why its 
contribution is so large in these massive, and prob- 
ably bulge-dominant galaxies. This is because mas- 
sive bulge-dominant galaxies at z — have many late- 
type progenitors, and that is more pronounced in less 
massive galaxies (Rank 3) than in Rank 1 galaxies. 
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Fig. 7. — Average mass evolution of stellar components in direct progenitors (upper) and the fraction of each component divided by 
the mean total stellar mass of direct progenitors at each redshift (bottom). The left, middle, and right panels show the mean evolutionary 
histories of the direct progenitors of Rank 1,2, and 3 galaxies, respectively. Each color code represents each stellar component as follows: 
(1) sky blue: merger accretion of stars formed in quiescent mode, (2) orange: in-situ quiescent star formation, (3) blue: in-situ starburst, 
and (4) red: merger accretion of stars formed in burst mode. 



This is a reflection of the progenitor bias discussed ear- 
lier fc.f.lFranx fe van DokkumlkOQlHGuo fe Whitd l2008l: 
iParrv et al.ll2009t: iKavirai et al.ll2009[) ~ 

Second, the stars formed in burst mode are an ex- 
treme minority (~ 2%) in these massive galaxies. This 
result is somewhat unexpected. In the hierarchical uni- 
verse, larger galaxies are generally a product of numerous 
mergers between galaxies. Mergers, especially between 
similar-mass galaxies (major mergers), cause a starburst, 
and so a naive expectation is that larger galaxies would 
have a large fraction of stars formed in burst mode. How- 
ever, it is not that simple. Starbursts are usually a re- 
sult of major mergers, and major mergers are extremely 
rare once galaxies become massive, that is, at low red- 
shifts. Major mergers between massive galaxies some- 
times occur at low rcdshifts, but they are most often dry, 
between early-type galaxies without much cold gas (c.f. 
van DokkunJ[200a: IBundv et al1l20"07l [200l IParrv et al.l 
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2009t) . Thus, it has been suggested that the fraction of 



stars formed in merger- driven burst mode may be less 
than a few per cent (e .g. I Aragon- Salamanca et al.lfl998l : 
IRodighiero etld1l2l)lll ). 

Third, merger accretion is the dominant channel in 
the most massive galaxies. Roughly 70% of the stel- 
lar components in the most massive group (bottom left 
panel) form outside direct progenitors and get accreted 
via mergers. This fraction is smaller for Rank 2 (40%) 
and Rank 3 (20%) galaxies, but still substantial. It is 
interesting to note that much of the mass difference be- 
tween ranks is attributed to the difference in the amount 
of merger accretion. For example, roughly 80% of the 
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Fig. 8. — Frequency of mergers (per Gyr per galaxy) with a mass 
ratio m2/mi > 0.1 experienced by direct progenitors. The solid, 
dashed, and dotted lines indicate the merger rates of Rank 1, 2, 
and 3 galaxies, respectively. 



mass difference between Ranks 1 and 2 comes from the 
differences in merger accretion. The value is ~ 60% be- 
tween Ranks 2 and 3. Massive galaxies are so mainly 
because they have a substantial amount of in-situ qui- 
escent star formation, but the most massive galaxies are 
so because they have acquired a large amount of mass 
through merger a ccretion. The result is s u pport ed by 
previous studies. I Aragon- Salamanca et al.l (119981 ) con- 
firmed that BCGs have experienced no or negative evo- 
lution in luminosity during < z < 1 from the K-band 
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Fig. 9. — The specific star formation rates (upper) and the spe- 
cific stellar accretion rates via mergers (bottom) of the direct pro- 
genitors in the three groups. See the text for their definitions. The 
solid, dashed, and dotted lines represent Rank 1, 2, and 3 galaxies. 

Hubble diagram for a sample of BCGs while they have 
increased their mass by a factor of two to four, depend- 
ing on cosmological parameters. Thus, it has been un- 
derstood that merger and accretion may be the most 
plausible expla nation for the evolut i on. U sing a semi- 
analytic model, iDe Lucia fc Blaizotl (|2007f ) showed that 
most stellar components in model BCGs are formed in 
the very early age of the Universe (80 % at z ~ 3) 
in small progenitors and accreted o nto BCGs f a r later 
(50% after z ~ 0.5) via mergers. IQser et al.l pOlOh 
presented similar results, using numerical simulations. 
About 80% of stellar components in simulated massive 
galaxies (M* > 2.4 x IQ^Mq at z=0) are formed outside 
in the early age (z > 3) and brought into the massive 
galaxies via mergers and accretion. The massive galax- 
ies double their mass after z ~ 1. They reve aled that the 
fractio n gets smaller in less massive galaxies. IParrv et al.l 
(|2009f ) also found a similar trend in their investigation on 
two separate semi-analytic models. They demonstrated 
that the contribution of mergers to the bulge growth ex- 
ceeds that of disk instability at M* > 10 115 M Q . Consid- 
ering the fact that such massive galaxies (M« > 10 11 M Q ) 
are likely early type (e.g. lBell et al. 2003), it implies that 
mergers play a more important role in the growth of mas- 
sive galaxies. It should however be noted that disk in- 
stability which is more effective to the smaller late-type 
galaxy evolution may play a role in such progenitor galax- 
ies of present-day massive early types. 
Figure [8] shows the merger rate evolution for bary- 
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Fig. 10. — Comparison of the mean specific star formation rates 
(solid lines) and the mean specific stellar accretion rates via merg- 
ers (dotted lines) of the direct progenitors of Rank 1 (upper), 2 
(middle), and 3 (bottom) galaxies. 

onic mass ratios greater than or equal to 1:10. Although 
merger rates show stochastic effects, there is a clear de- 
creasing tendency with time, whereas the star formation 
rates of the three groups drop more sharply, as shown in 
Figure [51 In general, more massive galaxies are likely to 
be involved in galaxy mergers more frequently, so that 
more massive galaxies have many more stellar compo- 
nents born outside and accreted via mergers, as discussed 
above. During the whole calculation, Rank 1 galaxies un- 
dergo about 9.0 mergers with a mass ratio greater than or 
equal to 1:10 while those in Ranks 2 and 3 experience 4.0 
and 1.0 mergers, respectively. During < z < 1, Rank 
1, 2, and 3 galaxies experience 3.5, 1.8, or 0.6 mergers for 
the same mass ratio criterion. This explains the higher 
contribution of merger accretion in more massive galax- 
ies, as illustrated in Figure [71 

4.3. Specific Star Formation Rates and Merger 
Accretion Rates 

We found in the previous section that in-situ star for- 
mation and merger accretion were the two most signif- 
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icant channels for the stellar mass growth of massive 
galaxies. In this section, we scrutinize their time evo- 
lution in greater detail. Specific star formation rates 
(SSFRs) are normalized growth rates of star formation 
histories. Likewise, we hereby define the "specific stellar 
accretion rate" (SSAR) to evaluate a normalized growth 
rate via mergers as follows: 



SSAR = 



AM 
M(t)At' 



(19) 



where M(t) is the mass of a galaxy at an epoch, and 
AM is an increment of mass by mergers during a time 
step At. Because we allow diffuse stellar components due 
to galaxy mergers, mass increment can be expressed as 
AM = (1 — jscatter)M<s Sa t, as described in Section 2.1. 

Figure M shows the evolution of the SSFRs (up- 
per) and of the SSARs (bottom) of direct progeni- 
tors. In general, more massive galaxies have lower 
SSFRs, as observations have shown (e.g. iSalim et al.l 
120071 iSchiminovich et al J 120071 ) . while they have higher 
SSARs than less massive galaxies. Star formation 
rates are decreased by the depression of gas cooling 
rates due to an increase in the cooling timescale via 
the g r owth of halos supernova feedback (IWhite fc Reesl 
[1971 IDekel fc Silld fl986t IWhite fc Frenkl 11991ft and/or 
AGN feedback dSilk fc Reesl I1998D . Furthermore, the 
cold gas reservoir of a galaxy could be reduced by feed- 
back processes. Besides, if a galaxy orbits around a more 
massive galaxy, it becomes redder as it loses its hot gas, 
which is a source of cold gas, and i ts cold gas reservoir 
by tidal and ram pressure stripping ( Gunn fc Gottl 1972j 
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20101). On the other 



hand, accretion of stellar components via mergers is de- 
termined by gravitational interactions between host and 
subhalos alone; hence, the accretion rate could remain 
relatively steady as halos continue to fall into other more 
massive halos over time. 

Figure [10] displays a comparison of the SSFRs and 
the SSARs of direct progenitors. As illustrated in Fig- 
ure quiescent star formation dominates the stellar 
mass growth history in (Rank 3) galaxies (bottom 
panel). In Rank 2 galaxies, they are comparable to each 
other most of the time. However, in the most massive 
(Rank 1) galaxies, merger accretion takes over star for- 
mation as the most important channel of stellar mass 
growth around z ~ 2 (top pa nel). Our r e sult i s qualita- 
tively consistent with that of lOser et al.l (|2010| ). 

5. SUMMARY AND DISCUSSION 

We have investigated the assembly history of stellar 
components of massive (super-L*) galaxies, using semi- 
analytic approaches. Our major results can be summa- 
rized as follows. 

• More massive galaxies grow in number faster than 
less massive galaxies, as a result of the hierarchical 
nature of galaxy clustering. This result is consis- 
tent with the recent observation of MK10. 

• The conflict between the predictions from hierar- 
chical models and downsizing is reconcilable. If we 
consider only direct progenitors of massive galaxies, 



our models suggest "upsizing" rather than downsiz- 
ing; that is, the direct progenitors of more massive 
galaxies grow more slowly. However, if we consider 
all the progenitors, direct and collateral, the com- 
bined mass suggests downsizing. Our models sug- 
gest that more massive galaxies have older stellar 
ages but younger assembly ages. 

• Merger-induced "bursty" star formation is negligi- 
ble compared to quiescent disk-mode star forma- 
tion despite the fact that massive galaxies form 
through numerous mergers. This is because most 
of the gas-rich major mergers occur at high red- 
shifts when galaxies are small, and recent major 
mergers tend to be rare and "dry" . 

• Merger accretion is a growingly more important 
channel of stellar mass growth in more massive 
galaxies. It accounts for 70% of the final stellar 
mass in the most massive galaxies in our sample 
(logM/A/ Q = 11.5 — 12.0). It is merger accretion 
that causes much of the mass difference between 
massive galaxies. This implies that environments 
play a central role in the growth of massive galax- 
ies. 

• In the most massive galaxies, which are likely 
brightest cluster galaxies, merger accretion has re- 
mained the most important channel of stellar mass 
growth ever since z ~ 2. 

The origin of massive galaxies is a pivotal subject of 
cosmological paradigms and subsidiary galaxy formation 
theories. The simplest views based on some pieces of 
observations, such as downsizing, may favor simplistic 
formation scenarios, while dynamical models of the uni- 
verse based on the current cosmology and other statis- 
tical aspects of observations, such as galaxy luminosity 
functions, indicate the other direction. This conflict is at 
its maximum when it comes to the formation of massive 
galaxies. The goal of this study is to reconcile the various 
perspectives and to have an accurate understanding on 
their formation. In this study, we showed that galaxy 
models in the hierarchical paradigm provide explana- 
tions to seemingly contradicting empirical constraints: 
the faster growth of the nu mber of more massi ve galax- 
ies (MK10) and downsizing (jCowie et al.lll996l ). This is 
encouraging. 

Galaxy formation models, whether hydrodynamic or 
semi-analytic, are still incomplete in many aspects: it is 
often claimed that much of their incompleteness is caused 
by our limited knowledge of baryon physics. The suc- 
cess of understanding massive galaxy formation on the 
other hand seems to be more hinged upon our knowledge 
of large-scale clusterings, and thus dark matter physics. 
Massive galaxies achieve their grandeur through mergers; 
thus, only by a realistic consideration of large-scale clus- 
tering information is it possible to accurately reconstruct 
their formation history. 
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